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We calculate the static dielectric tensor of a periodic system of aligned anisotropic 
dielectric cylinders. Exact analytical formulas for the effective dielectric constants 
for the E- and H- eigenmodes are obtained for arbitrary 2D Bravais lattice and 
arbitrary cross-section of anisotropic cylinders. It is shown that depending on the 
symmetry of the unit cell photonic crystal of anisotropic cylinders behaves in the 
low-frequency limit like uniaxial or biaxial natural crystal. The developed theory of 
homogenization of anisotropic cylinders is applied for calculations of the dielectric 
properties of photonic crystals of carbon nanotubes. 

PACS numbers: 
42.70. Qs, 41.20.Jb, 42.25.Lc 

I. INTRODUCTION 

Periodic dielectric structures - photonic crystals (PC) - have found many useful techno- 
logical applications.— Progress in the development of optical devices that operate using the 
principles of the control of light as proposed by Yablonovich^ gave rise to theoretical studies 
of the general properties of the spectra of PC's.^ In particular, the region of low frequencies, 
where the ideas and methods of the theory of homogenization^ are applicable, has attracted 
a significant amount of attention in the last decade. In the low-frequency limit the light 
wave in a non-absorbable periodic medium has linear dispersion, uj oc k. This allows the 



replacement a real inhomogeneous medium by an effective homogeneous one with dielectric 
permittivity 

In the general case, this effective parameter depends on the direction of propagation k = 
k/fc and has tensor structure. The latter property is emphasized for 2D PC's, which are 
anisotropic uniaxial or biaxial crystals, depending on the symmetry of the unit cell.— Unlike 
this, 3D PC's may be isotropic. 6 

Optical anisotropy of the PC's studied in Refs. f 5 ' 6,7 ] is determined by the geometry of 
the unit cell only. The constituents themselves are considered to be isotropic dielectrics. 
This is not the case, for example, for a structure of aligned carbon nanotubes arranged 
periodically in the plane perpendicular to the tubes. Here anisotropy manifests itself at the 
"microscopic" level, since the nanotubes ("atoms" of the PC) possess a natural anisotropy. 
This anisotropy originates from the layered structure of the crystal of graphite, which has 
different dielectric constants along the c-axis and in the perpendicular plane. The static 
values of these dielectric constants are e\\ = 1.8225 and e±_ = 5.226. 15 The elongated topology 
and the natural anisotropy of graphite cause PC's of carbon nanotubes to exhibit large 
optical anisotropy— Three-dimensional PC's with anisotropic dielectric spheres have been 
studied in Refs.fS]. It was shown, that, depending on the symmetry of the unit cell, the 
anisotropy of the spheres may be favorable for either broadening or narrowing the band 
gaps. 

High anisotropy of 2D photonic crystals may find interesting applications in nanophoton- 
ics as it was recently proposed by Artigas and Torner.^1 Namely, the surface of an anisotropic 
2D photonic crystal supports propagation of a surface wave predicted by Dyakonov.— The 
surface mode does not radiate and is localized close to the surface due to the interference 
between the ordinary and extraordinary waves. In natural crystals, it can be hardly ob- 
served because of the low anisotropy. Since it is a surface wave with very low energy losses, 
the Dyakonov wave may replace surface plasmons in the near-field optics and integrated 
photonic circuits. 

An effective medium theory for the PC of anisotropic carbon nanotubes was proposed in 
Ref. £12]. This theory is based on the Maxwell- Garnet approximation.— and it leads to a 
simple analytical formula for the effective dielectric constant. The formula is valid at low 



filling fractions. Since the Maxwell- Garnett approximation does not discriminate between 
periodic and non-periodic arrangements of the nanotubes, the in-plane anisotropy of the 
effective dielectric constant is lost. In other words, e e // in Eq. $IJ is independent on k, i.e., 
the Maxwell- Garnet approximation always leads to an effective medium which is equivalent 
to a uniaxial crystal. The PC considered in Ref. p^] has a square unit cell that gives rise to 
in-plane isotropy of the effective dielectric tensor. 

In this paper we extend the results of the theory of homogenization 5 to the case of 
anisotropic dielectric cylinders. Exact analytical formulas are obtained for the princi- 
pal effective dielectric constants of a 2D PC with an arbitrary cross-sectional form of 
anisotropic cylinders, arbitrary Bravais lattice, and filling fraction. We compare our re- 
sults with the results obtained for the PC of carbon nanotubes in the modified Maxwell- 
Garnet approximation.— This comparison shows that the Maxwell-Garnett approximation 
gives correct results for a dilute system with filling fraction less then 5%. However, at higher 
filling fractions the Maxwell- Garnet approximation overestimates the values for the effective 
dielectric constant. For a close-packed structure the error is about 10%. We also consider 
a PC with a rectangular unit cell and calculate two different in-plane dielectric constants. 
In this case the corresponding effective medium is a biaxial crystal. It was argued that 
for the if -polarized mode the effective dielectric constant for hollow and solid cylinders are 
practically indistinguishable.— Using our approach, we study the effect of the internal cavity 
on the effective dielectric constant and show explicitly how the effective dielectric constant 
decreases with the internal radius. 

II. THE FOURIER EXPANSION METHOD IN THE LONG- WAVELENGTH 

LIMIT 

We consider a 2D periodic structure of dielectric cylinders with their axes parallel to 
z and whose cross section can have an arbitrary shape. The cylinders are imbedded in 
a dielectric matrix. A 2D PC supports propagation of two uncoupled modes with either 
.E-polarization (where the vector E is parallel to the cylinders), or if-polarization (in this 
case the vector H is parallel to the cylinders). The background material is an isotropic 
dielectric with permittivity E\, and the cylinders are rolled up from an anisotropic dielectric 
sheet characterized by a tensor i^ a \ For carbon nanotubes, this tensor has two different 



eigenvalues, and in cylindrical coordinates is represented by a diagonal matrix with elements 
Egg = s'fj = e± and e!$ = en. As a whole, the periodic inhomogeneous dielectric medium is 



characterized by the coordinate-dependent dielectric tensor, 



f e xx (r) E xy (r) ^ 
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Inside the cylinders this tensor coincides with i^ a > and outside the cylinders it reduces to a 
scalar, £^8^. 

The wave equations for the E- and if -polarized modes with frequency u have the following 
form: 

V 2 E = ^E zz (r)E, (3) 

a ( dH\ oo 2 TT n 

dz[ a *w j ) + # H=0 > x *> x * = x >v- ^ 

Here E = E(x,y) and H = H(x,y) are the amplitudes of the E and H monochromatic 
eigenmodes, respectively, and is a 2 x 2 Hermitian matrix with determinant 1: 

a = 1 | £ ™( r ) £ *y( r ) | 

£xx\ r ) e yy\ r ) ~ £ xy( r )£yx( r ) \ E yx (r) E yy (r) j 

The determinant in the right-hand side can be written as a product of two eigenvalues. 
Within the graphite wall of the cylinders this product is £||£j_. Outside the wall it is either 
(e^) 2 (for r being outside the cylinders) or 1 for the interior region of the cylinder, which 
is free from the dielectric material. 

Eq.® depends only on the zz component of the dielectric tensor, i.e., it is insensitive to 
the in-plane anisotropy. Therefore the effective dielectric constant in the long- wavelength 
limit is the same as for the parallel arrangement of isotropic cylinders (not necessarily 
periodic). It is given by a simple formula 

£ ifl = £ zzi (6) 

where 

A c J A c 

is the average over the area A c of the unit cell zz component of the tensor (j2J. For a binary 
composite e\^ = e zz = fe± + (1 — f)£b, where / is the filling fraction of the component a. 
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To obtain the long-wavelength limit for Eq. we apply the method of plane waves.— 

Using the Bloch theorem and the periodicity of the function a^-(r), we get the Fourier 
expansions, 

H(r) = exp (zk • r) £ M G ) ex P (^ G • r )> ( 8 ) 
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where the Fourier coefficients ciij(G) are given by 
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Here G gives the reciprocal-lattice vectors. Indices a, b, and c at the integrals label the 
domains of integration - within the graphite walls, within the dielectric matrix, and within 
the interior of the hollow cylinders respectively. Substituting Eqs. (jHJ) and Q into Eq. (J1J) 
we get a generalized eigenvalue problem in G-space, 

$>;( G - G')(k+ G) 4 (k + G')^k(G') = (u 2 /c 2 )h k (G), i,j = x,y. (11) 

G' 

Eq. (jllj) is an infinite set of homogeneous linear equations for the eigenfunctions hk(G). 
The dispersion relation u = u> n (k) (n = 1, 2, . . .) is obtained from the condition that this set 
has non-trivial solutions. 

The periodic medium behaves as a homogeneous one if the Bloch wave (JHJ) approaches 
a plane wave. This occurs if the Fourier coefficients with G 7^ in (JHj) vanish in the long- 
wavelength limit. To obtain the behavior of /ik(G) we substitute G = in both sides of Eq. 
(fTTjl. divide the both sides by hk(G = 0) and take the limit as k — > 

1 = 1 £ M~ G ') fei G'j hi(G') . (12) 

Here a^- = ajj(G = 0) is the bulk average of the matrix (JSJ) and /i k (G) = /ik(G)//ik(G = 0) 
is the normalized Fourier coefficient. In the long-wavelength limit the coefficients of /i£(G') 
in the right-hand side are inversely proportional to k. In order to make the sum finite, the 
amplitudes of non-zero harmonics, /i£(G') must approach zero linearly with k. Thus, in the 
long-wavelength limit the Bloch wave (jHJ) can be written as a linear expansion over powers 
of k: 



H(r) = exp (ik • r) 



h + ]T hi(G) exp (iG ■ r) 
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Since the sum over G vanishes linearly with k, Eq. (j!3|) shows that the medium becomes 
homogeneous, i.e. the solution of the wave equation (JH) approaches a plane wave with an 
amplitude ho = h^G = 0) when k —>■ 0. 

Now, to calculate the effective dielectric constant we develop a theory of perturbation 
with respect to a small parameter ka (a is the lattice period). In Eq. f)l lj) we keep the linear 
terms and obtain the following relation, 

a, J (G)(;,k J + £ a tJ (G - G')GiG; hl(G') = 0. (14) 

The quadratic approximation is given by Eq. (|T2|) . which gives another linear relation 
between the eigenvectors /i£(G). Note that this relation is obtained from the eigenvalue 
problem Eq. (jll)) for G = and the linear approximation Eq. ([Tl^ is obtained for G ^ 0. 
The linear relations, Eqs. (fT2|) and (fTfj) . are the homogenized equations for the Fourier 
components of the magnetic field. These equations are consistent, if the corresponding 
determinant vanishes: 
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Aa tj (G - G^GiG'j - a lj (G)a mn (~G') G, n 3 n m G' n ] = 0. (15) 



Here n = ~k/k is the unit vector in the direction of propagation and A = (dijUiUj — £eff)- 
Since Eqs. (1T2J) and (fTljl are homogenous with respect to k, the dispersion equation (fT5|) 
depends only on the inverse effective dielectric constant, (u/ck) 2 . This fact is a manifestation 
of a general property: At low frequencies an electromagnetic wave has a linear dispersion in 
a periodic dielectric medium. 

Although Eq. ()15j) is an infinite-order polynomial equation in A, it turns out that it has 
only a unique nonzero solution. The fact that the second term in the determinant Eq. ()15j) 
is a product of two factors, one of which depends only on G and the other only on G', plays 
a crucial role. Omitting the mathematical details, which can be found in Ref . , we obtain 
the final answer for the inverse effective dielectric constant obtained from Eq. (|T3jl as: 
1 
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£ eff\ n > G,GV0 

Here [. . implies the inverse matrix in G-space. Eq. (jlfi)) is valid for an arbitrary form 
of the unit cell, geometry of the cylindrical inclusions, material composition of the photonic 
crystal, and the direction of propagation in the plane of periodicity. In the case when 
<2jj(G) = 77(G) 5ij Eq. (fTB^l is reduced to the formula obtained for isotropic cylinders.^ 



III. INDEX ELLIPSOID 



As any natural crystal, artificial PC in the long-wavelength limit can be characterized 
by an index ellipsoid.— Taking into account Eq. (JHJ) the equation for this ellipsoid can be 
written as follows: 

^I + M + io_ = i. (17) 

Here Xq, y$, Zq are three mutually orthogonal directions along which the vectors of the electric 
field, E, and of the displacement, D, are parallel to each other. For the E-mode we have 
E || D || z, i.e. the z direction coincides with z-axis. In the x — y plane the cross section of 
the index ellipsoid is given by Eq. (jl6|) . which can be rewritten in the canonical form as 

1 
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Here 

A ij = a ik( G ) a^-G^GkGl [a mn (G' - G)G m G' n ]~ l , i,j,k,l,m,n = x,y. (19) 

Eq. (fTSj) describes a rotated ellipse in the polar coordinates (p,f). The radius p{tp) = 
J e^J(<f) gives the index of refraction of H-mode and the angle ip is related to the direction 
of propagation, n = (co8(p,s\n.<p). The directions xq and yo coincide with the semi-axes of 
the ellipse given by Eq. (fTR|) and the in-plane indices of refraction y/e 1 , \fe 2 are given by 
the lengths of the semi-axes. The angle of rotation of the axes of the ellipse Eq. ()18|) with 
respect to the arbitrary axes x, y, is given by the relation 

9 A 

tan(20) = /y . (20) 
yy s*xx 

If the unit cell possesses a third- or higher-order rotational axis z, then the tensor is 
reduced to a scalar, A^ = A5ij, and the ellipse given by Eq. |T8|) is reduced to a circle. In 
this case the PC in the long-wavelength limit behaves like a uniaxial crystal; otherwise, it 
is biaxial. 



IV. UNIAXIAL AND BIAXIAL PC'S OF SOLID GRAPHITE CYLINDERS 



In this section, we study 2D PC of solid carbon cylinders arranged in square and rectan- 
gular lattices. In Cartesian coordinates the dielectric function of a carbon cylinder is given 
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by the tensor— 

/ Se|| + i5e± ^f(e||-e±) N 

Pf(e||-e±) Sqi + S^ • ( 21 ) 
V e ± j 

In numerical calculations we use Eq. (|T8|) and (|19|l . For rectangular and square lattices with 

circular cylinders the semi-axes of the index ellipsoid are directed along the basic lattice 

vectors. Because of the high symmetry of the unit cell, the off-diagonal elements of the 

tensor a^(G) vanish. The diagonal elements for hollow cylinders with outer and inner radii 

R and jR, respectively, (0 < 7 < 1) have the following form: 

' [I OT 1 + C 1 - 7 2 )] + H\l - - 7 2 )), G 
\gR (el 1 - e,- 1 ) [^(Gi?) - 1M1GR)} 

. + ( £ f - ^I 1 ) [ J o(Gi?) - J (7Gi?)]} , G ^ . 

The diagonal element a yy (G) is obtained from Eq. (j22j) by the replacement £j_ <->• en. 

In this section we consider solid cylinders, i.e. 7 = 0. The circles in Fig. ^ show the 

effective dielectric constant, given by Eq. (|T8|) . of the H-mode as a function of the filling 

fraction, / = ttR 2 /A c , for uniaxial PC with a square lattice. The number of G values (plane 

waves) considered in this calculation was 1200, which provided a good convergence in Eq. 

()18|) . The dielectric constant for the extra-ordinary mode (i?-mode), Eq. © (shown by 

triangles in Fig. ^) is always larger than that for the ordinary wave (if- mode). Therefore, 

the effective medium is a uniaxial positive optically anisotropic crystal. To check the validity 

of the Maxwell-Garnett approximation, we plot in Fig. ^ (squares) the effective dielectric 

constant proposed in Ref. p^] 
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(23) 



Here A = yje\\/e_\_. One can see that for all filling fractions the Maxwell-Garnett approx- 
imation gives overestimated values for the effective dielectric constant. For a very dilute 
system, / < 0.07, the Maxwell-Garnett approximation gives results that are practically in- 
distinguishable from the exact ones (See insert in Fig. QJ). For the close-packed array of 
cylinders the Maxwell-Garnet approximation overestimates the dielectric constant by about 
25%. 

In Fig. El we plot two principal dielectric constants for the biaxial PC of solid carbon 
cylinders with a rectangular unit cell. The ratio of the sides of the rectangle is 1:2. The 
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FIG. 1: In-plane effective dielectric constant for the H-mode for uniaxial PC of solid graphite 
cylinders with ey = 1.8225 and e± = 5.226 in air, = 1 (circles). Straight line (triangles) is 
the effective dielectric constant for the E'-mode, Eq. 0. The squares show the results of the 
Maxwell-Garnett approximation (|2rSJ). Insert shows the region of small filling fractions. 



difference between the two dielectric constants increases with the filling fraction, giving 
rise to a higher anisotropy of the corresponding effective medium. The Maxwell-Garnett 
approximation Eq. (|23|h which does not take into account the anisotropy of the unit cell, 
gives the values for emg that lie between the two principal values, S\ < Emg < £2- 



V. UNIAXIAL PHOTONIC CRYSTAL OF CARBON NANOTUBES 

In our model we consider the carbon nanotubes as hollow graphite cylinders. In the 
experimental study^ of the dielectric properties of carbon nanotues the outer radius of the 
cylinders was approximately R = 5 nm. The nanotubes formed a thin film and they were 
oriented along a specific direction. Although the nanotubes were not necessarily arranged 
periodically, one can assume that they formed almost a regular lattice, since the nanotube 
density is about 0.6 - 0.7 which is near the value of f c = 7r/4 w 0.785 for a close-packed 
structure. Thus, the separation between the nanotubes (the period of the square lattice 
d) slightly exceeds 2R, and in Ref. [12] it was estimated to be d = 10.15 nm. The inner 
radius 'jR = 0.25 — 2 nm.— was evaluated from the amount of electromagnetic absorption 
for the _E-polarized light. The four parameters /, R, 7, and A c = a 2 are not independent 
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FIG. 2: A plot of the principal effective dielectric constants for the PC of solid carbon cylinders 
arranged in a rectangular lattice. The larger (smaller) dielectric constant e\ (£2) corresponds to 
the direction of the vector E along the short (long) side of the rectangle. The Maxwell-Garnett 
dielectric constant is shown by the squares. 

but related by the formula, 



Substituting the aforementioned parameters of the square unit cell into this formula allows 
one to check that they are self-consistent. It is worthwhile to mention that the background 
material in the experiment^ is not air but the host material Delrin or Teflon with e& > 1. 
Since neither the density of the host material nor its dielectric constant is known, one cannot 
expect very good agreement between the experimental results^ and theory. In all theoretical 
considerations it was assumed that — 1. Because of this lack of experimental data, the 
effective medium theorie a 12 i 16 and the results shown in Fig. [I] give lower values for e e ff than 
that observed in the experiment.— 

It is obvious that the inner cavity reduces the permittivity of an isolated nanotube as 
compared to a solid graphite cylinder of the same size. It was argued^ that for a periodic 
arrangement the effect of the inner cavity is less than that for a single cylinder and even 
can be ignored, if the ratio between the inner and outer radii 7 does not exceed 0.4. This 
conclusion was supported by comparing the results of the Maxwell- Garnet approximation 
Eq. (j2Hj) and numerical band structure calculations. In Fig. El we plot the dielectric constant 
for a square lattice of hollow carbon nanotubes and compare the exact results obtained from 



f = nR\l- 1 2 )/A c . 



(24) 




0.20 0.25 0.30 0.35 0.40 0.45 

R/a 

FIG. 3: The plot of the effective dielectric constant for square lattice of carbon nanotubes versus 
the outer radius for tubes with different ratios of the inner and outer radii, 7 = 0.1,0.3,0.5,0.7. 
The exact results are shown by circles and the results of the Maxwell-Garnett approximation are 
shown by squares. 



Eqs. ()18p . (|19|h and (J221) (shown by the circles) with the results given by the Maxwell- 
Garnett approximation (squares). One can see that, for the same outer radius, the effective 
dielectric constant drops with an increase of the inner radius. Thus, if the outer radius is 
fixed, the dependence on the inner radius cannot be ignored, even in the Maxwell-Garnett 
approximation. However, the effective dielectric constants exhibits much less sensitivity 
to the internal radius if it is plotted as a function of the filling fraction, Fig. |3J In the 
Maxwell-Garnet approximation there is no dependence on the parameter 7, therefore, 
this approximation is represented by a single curve in Fig. HJ Here, only the total amount 
of the dielectric material is important, but not the topology of the cylinders. In the exact 
theory the effective dielectric constant depends on the details of the microstructure of the 
photonic crystal, but as far as the filling fraction is concerned, the topology plays a much 
less important role. Since the cylinder is uniquely determined by either two parameters out 
of three, R, 7, and /, the curves in Fig. 0]may cross each other. This means that at the 
crossing point the values of / and 7 correspond to the same hollow cylinder. This can be 
easily seen from Eq. (J24|) . 
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FIG. 4: A graph of the effective dielectric constant for a square lattice of carbon nanotubes versus 
filling fraction for tubes with different ratios of the inner and outer radii, 7 = 0.1, 0.3, 0.5, 0.7. The 
exact results are shown by circles and the Maxwell-Garnett approximation is shown by squares. 



VI. CONCLUSIONS 



We calculated the low-frequency dielectric tensor for 2D photonic crystal of anisotropic 
parallel cylinders arranged in a periodic lattice in the perpendicular plane. The exact ana- 
lytical formula for the principal values of the dielectric tensor was obtained. The results are 
applied for the periodic arrangement of carbon nanotubes which are rolled up from uniaxial 
graphite crystal with static values of the dielectric tensor e\\ = 1.8225 and e± = 5.226. It was 
shown that the interior (vacuum) region of the nanotubes has a small effect on the dielectric 
properties of the photonic crystal and can be ignored. Although we are interested in the 
static dielectric tensor, it is clear that the developed long-wavelength limit approach remains 
valid, even for optical frequencies since the period of the lattice of carbon nanotubes d — 10 
nm is much less than the optical wavelength A ~ 500 nm. To calculate the dynamic dielectric 
tensor, one has to substitute in the general formula Eq. (fbojl the corresponding frequency- 
dependent values for en and e±. Of course at finite frequencies Eq. (fTo|) gives the real part 
of the dielectric function. Calculations of the imaginary part require a generalization of the 
presented theory. This result will be reported elsewhere. 

The exact theory presented here allows a calculation of the effective dielectric constant 



of carbon nanotubes imbedded in a gas. Due to high absorbability of nanotubes, the con- 
centration of gas in the interior region of the nanotubes may be different from that in the 
atmosphere. This leads to slightly different dielectric constants of the material in the interior 
and exterior regions of the cylinders. This effect can be registered by precise measurements 
of the shift of the resonant frequency of a resonant cavity.-^ 8 * Thus, the proposed theory may 
find applications in the microwave detection of poisson gases in the atmosphere. 

One more interesting application of carbon nantonube photonic crystal is related to its 
huge anisotropy of the effective dielectric constant. Recently Artigas and Torner— demon- 
strated that the electromagnetic surface wave (Dyakonov waveii) can propagate along the 
surface of a photonic crystal with high optical anisotropy. This wave propagates in a loss- 
less dielectric medium and decays much slower than the surface plasmon. Since crystals 
with huge optical anisotropy are rare in nature, carbon nanotube photonic crystals may be 
considered as a promising material for integrated photonic circuits. 
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